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ABSTRACT 

We examine a sample of 48 Swift/WOT long Gamma-ray Burst light curves and find 
a correlation between the logarithmic luminosity at 200s and average decay rate deter- 
mined from 200s onwards, with a Spearman rank coefficient of -0.58 at a significance 
of 99.998% (4.2a). We discuss the causes of the log L 2 oo s — «>200s correlation, find- 
ing it to be an intrinsic property of long GRBs, and not resulting from the selection 
criteria. We find two ways to produce the correlation. One possibility is that there 
is some property of the central engine, outflow or external medium that affects the 
rate of energy release so that the bright afterglows release their energy more quickly 
and decay faster than the fainter afterglows. Alternatively, the correlation may be 
produced by variation of the observers viewing angle, with observers at large viewing 
angles observing fainter and slower decaying light curves. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are intense flashes of gamma- 
rays that are usually accompanied by an afterglow, longer 
lived emission that may be detected at X-ray to radio wave- 
lengths. Our understanding of GRB X-ray and optical/UV 
afterglows was revolutionised by the launch of Swift, a satel- 
lite dedicated to the detection of GRBs and observation of 
their 7-ray, X-ray and optical/UV emission (Gehre ls et alj 
120041 ). The best studied subclass of GRB are the long GRBs 
(LGRBs), which h ave observed prompt 7-r ay emission du- 
rations of > 2s l|Kouveliotou et all [1993). Several stud- 
ies have investigated the X-ray emission of LGRBs, using 
large samples, to id entify characteristic t emporal and spec- 
tral behaviours (i.e iNousek et~aH l200d : lEvans et all 120071 
2009). Similar investiga tions have been perfo r med at opti- 
cal/IR wavelengths (i. e iMelandri et al.l 120081 ; iRvkoff et al.1 
12009] ; "Oates et aU2009l h but these have tended to have much 
smaller samples, due to the lower detection rate of the op- 
tical emission in comparison to the X-rays (the detection 
rate for Swift's X-ray and ultr a-violet optical telescopes are 
~ 96 % (|Burrows et al.l 120081 ) and ~ 29% (|Roming et al.1 
2009), respectively for LGRBs). The low detection rate is 
generally att ributed to extincti on due to high levels of dust 
in the host (|Fvnbo et al.|[200lT ). and/or to high redshift at 
which the optical/UV emission will be absorbed by neu- 



tral hydrogen along the li ne of sight (|Groot et alJI 199*3 ) ; see 
also iGreiner et al.1 (|201ll ) for a recent study on the cause 
of optically dark GRBs. However, Swift has now observed 
over 600 GRBs and we are now in a position to collate a 
large number of well-sampled IR/optical/UV LGRB after- 
glows from which we can draw a representative picture of 
th eir behaviour and collective prop erties (see recent papers 
bv lKann et"ai1l2010l ; lLi et al.ll2012l ). 



In lOates et alj (|2009l ). we performed a statistical in- 
vestigation of 26 optical/UV LGRB afterglows, finding a 
correlation between the observed «-band magnitude at 400s 
and the average UVOT light curve decay rate determined 
from 500s. We also tested for an equivalent rest frame cor- 
relation, but, due to the small sample size, this could not 
be confirmed or excluded. Here we use a larger sample of 
48 high quality LGRB UVOT light curves to re-examine if 
there is a correlation between optical/UV afterglow intrinsic 
brightness and light curve decay rate. 

This paper is organized as follows. In § [5] we describe 
the data reduction and analysis. The main results are pre- 
sented in § [3] Discussion and conclusions follow in § [4] and 
[S] respectively. All uncertainties throughout this paper are 
quoted at la. The temporal and spectral indices, a and /3, 
are given by the expression F(t, v) oc t a v^ . 
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Figure 1. Optical luminosity light curves of 56 GRBs at rest- 
frame 1600A. For clarity, 3 a upper limits are not included. 
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2 DATA REDUCTION AND ANALYSIS 
2.1 SAMPLE 

Our sample began with 69 LGRB s from the second Sw i 
UVOT GRB afterglow catalogue (|Roming et al.l ITn prep 
which were observed between April 2005 and December 
2010. The sample was selected using the same criteria as 
lOates et al.l (|2009l ): the optical/UV light curves must have 
a peak UVOT u-band magnitude of ^ 17.89 (equivalent to 
1 ct s _1 ), UVOT must observe within the first 400s until 
at least 10 s after the BAT trigger and the colour of the 
afterglows must not evolve significantly with time, meaning 
that at no stage should the light curve from a single filter 
significantly deviate from any other filter light curve when 
normalized to the v filter. Only GRBs 060218 and 100814A 
were excluded as we considered them to have strong colour 
evolution. These criteria ensure that a high signal-to-noise 
(SN) light curve, covering both early and late times, can be 
constructed from the UVOT multi-filter observations. 

In order to obtain the best SN light curve for each GRB, 
a single light curve was constructed from the mu l ti-filte r 
light curves, following the method in lOates et alj {2009). 
The main steps are to normalize the multi-filter light curves 
to the v filter and then to gro up them toge t her u sing a bin- 
size of At/t = 0.2. Following fpates et al.l l|2009h . for each 
GRB, we take the onset of the prompt 7-ray emission as 
the start time of the UVOT light curve. Since the BAT may 
not trigger at the start of the prompt emission, we take the 
start time of the Tgn paramete o be the start time for the 
UVOT light curve. 



2.2 Luminosity Light curve 

Luminosity light curves were produced for all GRBs with 
spectroscopic or photometric redshifts in the literature and 
for which host E(B-V) values could be determined. For a 
further 3 GRBs we were able to derive photometric red- 
shifts from joint XRT-UVOT Spectral Energy Distributions 



1 The Tgo parameter is determined from the gamma-ray event 
data for each GRB, by the BAT processing script. The results 
of the processing are publicly available and are provided for each 
trigger at http://gcn.gsfc.nasa.gov/swift_gnd_ana.html 



(SEDs), following the method of lKriihler et all l|201lh : z = 
1.2 ±0.1 for GRB 060510A, z = 3.1 ±0.1 for GRB 090401B 
and z = 1.85 ± 0.10 for GRB 100805A. For 13 GRBs, no 
redshift was available nor could we derive a photometric 
redshift. 

For each of the 56 GRBs with host E(B-V) and redshift 
measurements, we converted the single filter count rate light 
curve to luminosity at a common rest frame wavelength. In 
order to select the common wavelength and determine the 
resulting k-correction factor for each light curve, an SED 
was computed for each GRB following the methodology in 
lOates et all {2009) an d using the count rate -to-flux conver- 
sion factors given in iBreeveld etaD (|201lh . The common 
rest frame wavelength was selected to maximise the number 
of GRBs with SEDs that covered this wavelength and to 
be relatively unaff ected by host extinction. As was found in 
lOates et al.l {2009T I. the wavelength that best satisfies these 
conditions is 1600A. For each GRB, the k-correction factor, 
k, was taken as the flux density at the wavelength that cor- 
responds to 1600 A in the rest frame, -F1600, divided by the 
flux density at the observed central wavelength of the v filter 
(5402 A), F v , which was multiplied by (1 ± z), where z is the 
redshift of the GRB such that k = (F 1600 /(F V * (1 + z))). 
For those GRBs with SEDs not covering 1600A, an average 
k value was determined from the other GRBs in the sample, 
which have SEDs covering both 1600A and the v filter rest 
frame wavelength. 

To convert to luminosity, the count rate light curves 
were initially corrected for Galactic extinction and then con- 
verted into flux per unit frequency. These flux densities were 
then converted to luminosity at 1600A using: 



L(1600) = AnD 2 L F v k 



(1) 



where L(1600) is the luminosity at a 1600A and Dl is the 
luminosity distance. Finally, the luminosity light curves were 
corrected for host extinction, w hich was calculated from 
the host A v values reported in Schady et all {2010). For 
those not reported in lSchadv et al] ~i 2010t) . we used the same 
method to determine the host A v . 



3 RESULTS 

In Figure [T] we show the luminosity light curves at 1600 A, 
in units of erg s _1 Hz -1 . The light curves are clustered in a 
single group, having the largest range in luminosity at the 
earliest epochs, which becomes narrower as the light curves 
decay. This characteristic is also seen in the light curves 
given m iKann et al.l (|2010T) . The narrowing of the range in 
luminosity with time suggests that the most luminous GRBs 
decay the quickest and the less luminous GRBs decay more 
slowly. In the next section, we determine if there is a sta- 
tistically significant correlation between the brightness and 
decay rate and quantify this correlation. 



3.1 Luminosity Decay Correlation 

In order to test for a correlation we determined the intrinsic 
brightness at 200s and the average decay rate determined 
from 200s onwards. To determine the intrinsic brightness, 
we used the IDL interpolation function, Interpol, on the 
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Figure 2. Average decay index determined from the luminosity 
light curves after 200s versus luminosity at 200s. The red solid 
line represents the best fit regression and the blue dashed line 
represents the 3<r deviation. 



data between 100 and 2000s and interpolated the logarith- 
mic luminosity at 200s, log L2oos- For the average decay 
rate, cv>200s, using error weighted least squares we fit a sin- 
gle power law to the light curves from 200s onwards. A rest- 
frame time of 200s was chosen as by this time all the light 
curves in the sample have observations. For six GRBs, there 
were too few data points between 100 and 2000s to constrain 
the luminosity at 200s and we excluded two further GRBs 
for which the error on the decay index was greater than 3 
times the sample mean decay error. The resulting values for 
log L200S and a>200s are shown in Fig. [2] A Spearman rank 
test of log L200S against q>200s, gives a correlation coefficient 
of -0.58 at a significance of 99.998% (4.2a). This indicates 
that log L200S and a>2oos are statistically correlated and 
confirms that luminous optical/UV afterglows decay more 
quickly than less luminous ones. 

In Fig. [3] we show the distributions of the logarithmic 
luminosity at 200s and 6 hours, log Lghrs- If the brighter op- 
tical/UV afterglows decay more quickly than the fainter ones 
then at late epochs the luminosity distribution must become 
narrower and the correlation should become weaker and/or 
insignificant. We observe both of these effects in our sample. 
In Fig. [3] the 6hr luminosity distribution has a standard de- 
viation of 0.69 which is smaller than the standard deviation 
of the 200s lumin osity distribution, which is 0.84 (see also 
iKann et alj|2010h . and, a Spearman rank test of log L6hrs 
and a > 200s indicates no significant correlation with a coeffi- 
cient of 0.01 and a significance of 8%. 

In order to quantify the relationship between log L200S 
and «>200s, we performed a linear regression using the IDL 
routine, fitexy, which takes into account the errors on both 
parameters. This analysis provides a linear relationship of 
log Laoos = (-3.636 ± 0.004)o? + (28.08 ± 0.13), given as a 
solid red line in the left panel of Fig. [5] All GRBs reside 
within 3x the rms deviation, indicated by blue dotted lines. 

We looked to see if the log L200S — a>200s correlation was 
due to both parameters being related to redshift, 2. A Spear- 
man rank correlation between these parameters gives: a co- 
efficient of 0.62 at a significance of 4.7a for log L200S — z, and 
a coefficient of -0.33 at a significance of 98% for a>200s — z. 
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Figure 3. Distribution of logarithmic luminosities at rest frame 
200s (top) and 6hrs (bottom). 



The correlation between log L200S and redshift is expected 
as we are able only to detect GRBs at high redshift if 
they are very luminous. The correlation between cv>200s 
and redshift is much weaker than the correlation between 
log L 2 oos and q>2oos suggesting that the L 20 o s - «>200s cor- 
relation is not due to the implicit correlation between these 
two parameters and redshift. To confirm this, we calculate 
the partial Spearman rank correlation, which measures the 
degree of correlation between two parameters, L200S and 
Q>2 oos, excluding the effect of a third, in this case redshift 
(see lKendall fc Stuartlll979l . for further details). Using this 
method we obtain a correlation coefficient of -0.50 with a 
confidence of 99.97% (3.5a), only a small reduction in the 
correlation coefficient from the standard Spearman rank cor- 
relation. This indicates that the correlation between L200S 
and a>2oos is not a result of the implicit correlation between 
these two parameters and redshift. 

The log L200S — a>200s correlation could be due to 
chance or be a result of selecting the sample by their ob- 
served frame properties, specifically the exclusion of LGRB 
optical afterglows that were fainter than 17.89 in UVOT 
n-band and/or were not observed within 400s after the trig- 
ger. In order to eliminate these possibilities, we performed a 
Monte Carlo simulation. For each of a total of 10 6 trials, we 
simulated a distribution of 48 pairs of L200S — Q>2oos data 
points selected at random from linear distributions of L200S 
and a>200s, which have the same ranges as the observed 
sample. For each pair of data points, we produced a syn- 
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Figure 4. Optical/UV temporal and spectral indices for the sam- 
ple of 48 GRBs. The lines represent 3 closure relations and a 
colour scale is used to display the range in luminosity at 200s, 
L200s- 

thetic observed frame light curve using L200S and Q>200s and 
randomized values for redshift, host and Galactic extinction 
and p- where, ft, wa s used to determine the k -correction fac- 
tor (k = (1 + z) 1+/3 ; lLamb fc Reichartll200d ). Similar to the 
simulated L200S and a>200s points, these parameters were 
sampled randomly from linear distributions, which have the 
same range as the observed light curves. In order to sim- 
ulate the time it takes for Swift to point its narrow field 
instruments at the GRB location, we selected at random an 
observed frame light curve from our sample and used the 
time sampling of this light curve, i.e photometry times and 
durations, as the time sampling for the simulated observed 
frame light curve. We then determined if the simulated ob- 
served frame light curve met our selection criteria. For those 
that did not meet the selection criteria, we discarded the 
corresponding L200S — cv>200s data point and drew a new 
pair of values until the selection criteria were met. Once 48 
pairs of L200S and a>200s had been found with synthetic ob- 
served frame light curves that met the selection criteria, a 
Spearman rank correlation was performed on the simulated 
L200S — «>2oos distribution. 

Of the 10 6 trials, only 34 have a correlation coefficient 
equal to or indicating a stronger correlation than the real 
L200S — cv>200s distribution. This indicates that, at 4.1a con- 
fidence, the L200S — Q>2oos correlation is not due to our se- 
lection criteria nor does it occur by chance and therefore 
implies that the L200S — cv>2oos correlation is intrinsic to 
LGRBs. 



4 DISCUSSION 

Using 48 optical/UV light curves, we have found a signifi- 
cant correlation between the intrinsic brightness and average 
decay rate of LGRB o ptical/UV afterglows. A simila r cor- 
relation was found by iPanaitescu fc Vestrandl {2008) , who 
found for those GRBs with an observed rise in the opti- 
cal light curves, there was a correlation between peak flux 
and post- peak power-law decay r ate. T he correlation re- 
ported by IPanaitescu fc Vestrandl l|200ct ) is limited to only 



those GRBs which are observed to rise. Since the major- 
ity o f rising afterglows cease rising by ~ 400s l|Oates et alj 
I2009T ) . which is typically earlier than restframe 200s (i.e 
~ 400s/(l + z) < 200s), the log L200S — ">200s correla- 
tion determined in this work can be applied to optical after- 
glows with all types of early behaviour. Furthermore, as the 
UVOT typically begins observing within the first 100s after 
the BAT trigger, the majority of GRBs observed by Swift, 
with a detected optical counterpart, can be included in the 
correlation reported in this paper. In the following, we shall 
examine the possible ways to produce a correlation between 
log L200S and cv>200s- 

4.1 Is the correlation predicted by the standard 
afterglow model? 

4-1.1 The basic model 

The log L200S — a>200s correlation may be a natural result 
of the jet interacting with the external medium, producing 
synchrotron emission. Here we assume an isotropic, colli- 
mated outflow which is not energy injected. In this model, 
the luminosity is related to ct and ft by L v cx F v cc t a v^ , 
and a and are related linearly by the closure relations 
fe.g ISari. Piran fc Naravanl Il998h . defined by the external 
medium density profile and the ordering of the synchrotron 
frequencies. For the GRBs in this sample, due to our selec- 
tion criteria requiring no colour evolution, we may conclude 
that /3 does not vary during the course of our observations. 
We now examine two scenarios to determine if the correla- 
tion is a result of the basic afterglow model. 

In the simplest scenario, all optical afterglows arise from 
a single closure relation, in which a is not a fixed value. 
In this scenario, a and ft are related linearly and therefore 
a correlation between log L200S and /? should be expected. 
However, only a weak correlation is observed between /3 and 
q and there is no evidence for a correlation between /3 and 
log L200S (see Fig[4]), with Spearman rank coefficients of 0.26 
(92%) and -0.15 (68%), respectively. This scenario is there- 
fore not likely the cause of the log L200S — cv>200s correlation. 

However, we do not expect all optical afterglows to be 
on the same phase of the synchrotron spectrum. In the 
second scenario, we assume that the log L200S — cv>200s 
correlation is a result of more than one closure relation, 
which requires the use of multiple spectral segments and 
different density profiles. The spectral segments most com- 
monly found to satisfy optical afterglow production are: 
v m < v op t < v c and v c < u op t, where v c is the cooling fre- 
quency and u m is the peak frequency. The v m < v ovt < v c 
spectral segment requires different closure relations for when 
the ejecta interacts with a constant density medium and a 
wind- like density medium. The spectral segment, v c < u op t, 
is independent of the external medium density profile. The 
expected relations between a and /3 for these scenarios are 
indicated by lines in Fig. [4] If the log L200S — a>2oos cor- 
relation was produced by the optical afterglows resulting 
from multiple closure relations, in Fig. U we would observe 
the a and /3 data points with similar luminosities to cluster 
around a given closure relation. In Fig. [4] the data points 
do not appear to display this trend and therefore we find 
it unlikely that the basic standard model is causing the 
log L200S — Q>200s correlation. 
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4-1-2 Complex afterglow model 

The afterglow model is likely to be more complex than we 
previously assumed and there may be some mechanism or 
parameter that regulates the energy release in GRB after- 
glows and their decay rate. To satisfy our observations, this 
must occur in such a way that when the energy is released 
quickly the result is an initially bright afterglow which de- 
cays quickly. Conversely, if the energy is released slowly over 
a longer period, the afterglow will be less bright initially and 
decay at a slower rate. This may indicate that there is a nar- 
row range of energy provided to the outflow. One possible 
way to regulate energy release could be continued energy in- 
jection. If the central engine does not initially release all its 
energy, but releases it over a much longer period, the result 
could be a fainter afterglow which decays slowly. 

4.2 Observing angle, jet structure and degree of 
collimation 

The L200S — «>200s correlation may instead be due to a 
range in observing angle (i.e observer's angle relative to the 
jet axis), 9 a bs, with the fainter optical afterglows being ob- 
served at larger observing angles ( e.g Ramirez-Ruiz et alj 
l2005l;|Panaitescu fc Vestrandll2008h . |Panaitescu fe Vestrandl 
1 2008) show in their fig. 3. that, for a jet with uniform ve- 
locity distribution, an off-axis observer, i.e 9 bs/9jet > 1, 
where 9j et is the jet opening angle, will observe a shallower 
decay and observe the afterglow to be less luminous in com- 
parison to an observer who is observing closer to the edge 
of the jet, i.e o bs/9jet ~ 1- This effect should be observed 
for both constant density and wind-like external media. This 
model is more complex if structured outflows are considered. 
For structured outflows, off-axis viewers will also observe a 
shallower and fainter light curve in comparison with on-axis 
observers, but the convergence time and the range of decay 
rates will vary, depending on how the outflow is structured. 

A couple of tests may provide support for the log L200S — 
Q>200s correlation resulting from jet structure and observer 
viewing angle. First, we should expect to see convergence 
of the light curves at late times to a similar decay rate 
for all observing angles. And second, since afterglows that 
are viewed more off-axis will rise later, we should also ob- 
serve a correlation between afterglow brightness and peak 
time, though this may be c omplicated by jet structure. 
iPanaitescu fc Vestrandl l|2008h tested for this correlation us- 
ing 11 optical light curves with rises and find a strong cor- 
relation between peak luminosity and peak time consistent 
with this hypothesis though their test c ould not be applie d 
to the GRBs without observed rises fcf lKann et"aLll2010h . 

5 CONCLUSIONS 

We computed luminosity light curves at 1600A for 48 opti- 
cal/UV GRB afterglows. We find a correlation between lu- 
minosity at 200s and average decay rate from 200s onwards 
with a significance of 99.998% (4.2a). Regression analysis 
indicates a linear relationship between decay rate and lumi- 
nosity of log L 2 o 0s = (-3.636 ± 0.004)a + (28.08 ± 0.13). 

We used a Monte Carlo simulation to determine, at 4.1a 
confidence, that the L200S— a>200s correlation is intrinsic and 
not due to chance or our selection criteria. We determined 



that this correlation is not likely to be a natural consequence 
of the basic synchrotron afterglow model. Instead we find 
two possible ways to produce the correlation. The first is 
that there is some property of the central engine, outflow 
or external medium that affects the rate of energy release 
and rate of light curve decay, in such a way that for brighter 
afterglows the energy is released more quickly and decays 
more rapidly than the fainter afterglows. Alternatively, the 
correlation may be produced by a range in observing angles, 
with observers at large viewing angles witnessing fainter and 
slower decaying light curves. Understanding the origin of 
this correlation will have important consequences on our un- 
derstanding of the physics and geometry behind GRBs. 
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